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Abstract 

Hepp and Speer sectors were successfully used in the sixties and seven- 
ties for proving mathematical theorems on analytically or /and dimen- 
sionally regularized and renormalized Feynman integrals at Euclidean 
external momenta. We describe them within recently developed strate- 
gies of introducing iterative sector decompositions. We show that Speer 
sectors are reproduced within one of the existing strategies. 
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1 Introduction 



The so-called alpha representation [1-11] was initially used to introduce dimensional 
regularization [121 [13] and to prove various mathematical results on Feynman integrals 
O m El m [H [9l [ini E] . The standard way to analyze convergence of Feynman integrals 
is to decompose the initial integration domain of alpha parameters into appropriate 
subdomains (sectors) and introduce, in each sector, new variables in such a way that 
the integrand factorizes, i.e. becomes equal to a monomial in new variables times 
a non-singular function. This procedure turned out to be successful for Euclidean 
external momenta, i.e with g^)^ < for any partial sum, when the sectors of Hepp 
[3] and Speer [5] were introduced. 

However, in practice, one often deals with Feynman integrals on a mass shell 
or at a threshold. In this case, Hepp or Speer sectors generally do not provide a 
factorization of the integrand so that the analysis of convergence fails within this 
technique. Therefore, general theorems on such 'physical' Feynman integrals could 
not be proved up to now. 

Recently Binoth and Heinrich introduced sector decompositions of a new kind 
[Ti] . (See [15] for a review.) They provided a powerful method of evaluating Feynman 
integrals numerically in situations with severe UV, IR and collinear divergences. In 
contrast to Hepp and Speer sectors, the sectors of [14] are introduced iteratively, 
according to so-called sector decomposition strategies. The corresponding algorithm 
was implemented on a computer. Although this algorithm was successfully applied to 
numerically evaluate complicated Feynman integrals and to check existing analytical 
results (see, e.g. [HI [T71 [18]) it was unclear where this iterative procedure stops at 
some point, i.e. results in the factorization of the integrands so that one can apply it 
for numerical evaluation. Indeed, in some examples, closed loops appear within this 
algorithm. 

The first algorithm guaranteed to terminate was developed by Bogner and Weinzierl 
[in]. More precisely, certain strategies within this algorithm guarantee that closed 
loops do not appear. The algorithm works at least if squares (^Q'j)^ of partial sums 
of the external momenta are either negative or zerojj (We will formulate a more 
general condition in the next section.) The corresponding computer code is public. 
The second public algorithm [20] called FIESTA provided one more strategy (Strat- 
egy S) for sector decompositions which also leads to a factorization of the integrand 
for general Feynman integrals. It was successfully applied in [21|. 

The purpose of this paper is to describe Hepp and Speer sectors in an iterative 
way, within the modern technique of sector decompositions. Using examples of three- 
and four-loop diagrams, we will also see that Speer sectors are reproduced within 
Strategy S. In the next section, we describe the alpha representation and Hepp and 
Speer sectors in a simplified form. Then, in Section 3, we establish their connection 

^ Let us stress that Hepp and Speer sectors generally do not provide a resolution of the singular- 
ities in the parameter of dimensional regularization if some sum ^ qi is light-like. 
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with iterative sector decomposition. Finally, we discuss some perspectives. 



2 Hepp and Speer sectors 

For a Feynman integral with standard propagators (of the l/{rn? — k"^ — iO)"' form) 
corresponding to a connected graph F, the alpha representation has the following 
form: 

j^a+h{l-d/2)^hd/2 

Fr(gi, . . . ,g„; rf; ai . . . jOl) = 



X 



/OO /"OO 
... j J] a;'-iW-'^/2e^^/"-*S™?-'d«i . . . d«i 



where L and h is, respectively, the number of lines (edges) and loops (independent 
circuits) of the graph, n + 1 is the number of external vertices, a = ^ a^, and 

u = En«^' (2) 

In ([2]), the sum runs over trees of the given graph, and, in (JSj), over 2-trees, i.e. maxi- 
mal subgraphs that do not involve loops and consist of two connectivity components; 

is the sum of the external momenta that flow into one of the connectivity com- 
ponents of the 2-tree T. (It does not matter which component is taken because of 
the conservation law for the external momenta.) The products of the alpha param- 
eters involved are taken over the lines that do not belong to the given (2-)tree T. 
The functions U and V are homogeneous functions of the alpha parameters with the 
homogeneity degrees h and h + 1, respectively. 

An overall integration can be performed to obtain another well-known parametric 
represent atioE0: 

U7T^/^)''T{a-hd/2) 
Fr{qi,...,qn]d]ai...,aL) = — 



X 



n.r(aO 

^a-(h+l)d/2 



According to the well-known folklore Cheng- Wu theorem one can choose any sum of 
the alpha variables in the argument of the delta function. 



^ So, the code of [T^ works at least if each monomial of the a-variables in the function —V + 
U ^ mfai enters with a negative coefhcient. This condition is a little bit relaxed within the code of 
[20j where combinations of the type (a^ — aj)'^ are also admissible. 
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We imply that the graph F is a connected graph, i.e. any two vertices of T can be 
connected by a path in F. However, we are going to consider various subgraphs of the 
graph and they can be disconnected, i.e. consist of several connectivity components. 
A subgraph 7 of F is determined by a subset of lines and includes all the vertices 
incident to these lines. (Sometimes isolated vertices are added to a subgraph. For 
example, Mathematica produces isolated vertices as bi-connected components.) The 
number of loops of a subgraph equals to 

M7)=L(7)-V(7) + c(7), 

where V{'y) and 0(7) are, respectively the numbers of the vertices and connectivity 
components. 

An articulation vertex of a graph F is a vertex whose deletion disconnects F. Any 
graph with no articulation vertices is said to be biconnected (or, one-vertex-irreducible 
(IVI)). Otherwise, it is called one-vertex-reducible (IVR). In other words, in a IVR 
graph, one can distinguish two subsets of its lines and a vertex (an articulation vertex) 
such that any path between vertices from these two subsets goes through this vertex. 
From now on let us suppose that we are dealing with a IVI graph. It is natural to 
treat a single line as a IVI graph since we cannot decompose it into two parts. 

Any subgraph can be represented as the union of its IVI components, i.e. maximal 
IVI subgraphs. Consider, for example, the two-loop self-energy diagram of Fig. [H 
The subgraphs {1, 2, 5} and {1, 2, 3, 4} are IVI. The subgraph {1, 2, 3, 5} is IVR and 




Figure 1: A two- loop propagator diagram 

its IVI components are {1,2,5} and {3}. The subgraph {1,2,3} is IVR and its IVI 
components are {1}, {2} and {3}. 

A set / of IVI subgraphs is called an ultraviolet (UV) forest if the following 
conditions hold: 

(z) for any pair 7, 7' G /, we have either 7 C 7', 7' C 7 or £(7 fl 7') = 0; 
(a) if 7^, . . . , 7" G / and £(7*n7-') = for any pair from this family, the subgraph 
UiY is IVR. 

In other words, the number of loops in UiY (where 7* are disjoint with respect 
to lines and belong to a UV forest) is equal to the sum of the numbers of loops of 
7*. The term "UV" is used because the UV divergences are due to the integration 
over small values of ai where the exponent in ([T]) is irrelevant and they are generated 
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by the singularities of the factor U^^^'^. We are going to show that the resolution 
of the UV singularities can be performed by the use of sectors associated with IVI 
subgraphs. 

For example, the set {1}, {2}, {3} of subgraphs of Fig. [1] is a UV forest and 
{1,2,5}, {3} is also a UV forest but the set {!}, {2}, {3}, {4} is not a UV forest 
because the condition (ii) breaks down. 

Let be a maximal UV forest (i.e. there are no UV forests that include JF) of a 
given graph F. An element 7 G JF is called trivial if it consists of a single line and is 
not a loop line. Any maximal UV forest has h non-trivial and L — h trivial elements. 

Let us define the mapping a : JF — > £ such that (7(7) G >C(7) and (7(7) ^ -^^(7') 
for any 7' C 7, 7' G JF. The inverse mapping : C ^ T exists and can be defined 
as follows: o'^iV) is the minimal element of the UV forest T that contains the line I. 
Let us denote by 7_|_ the minimal element of T that strictly includes the given element 
7- 

For a given maximal UV forest JF, let us define the corresponding sector (/-sector) 

as 

Vjr = {a\ai < aa{j), / G 7 G JF} . (5) 

The intersection of two different /-sectors is of measure zero; the union of all the 
sectors gives the whole integration domain of the alpha parameters. For a given 
/-sector, let us introduce new variables labelled by the elements of JF, 

ai= n ' 

where the corresponding Jacobian is ^7 '^^^ inverse formula is 

^ _ f a;o-(7)/«o-(7+) if 7 is not maximal 
^ \ Q;o-(7) if 7 is maximal 

Consider, for example, the following maximal UV forest T of Fig. [1] consisting 
of 7I = {1}, 72 = {2}, 73 = {3}, 7^ = {1,2,5}, 7^ = F. The mapping a is 
cr(7^) = 1, o"(72) = 2, (7(73) = 3, o"(7'^) = 5, cr(7^) = 4. The sector associated with 
this maximal UV forest is given by Djr = {ai,2 < «5 < «4, «3 < 0^4} and the sector 
variables are t^i = ai/a^, t^2 = 02/0(5, t^s = a^/a^, t^i = a^/a^, t^s = a^. 

All the maximal UV forests of the given graph can be constructed at least in two 
ways. 

Way 1. We imply that the lines are enumerated. Let us consider the sequence 
of subgraphs 7^ consisting of lines {1, 2, . . . , /}, respectively, with / = 1, . . . , L. For 
each Z, let us take the IVI component of 7^ that includes the line I. The set of all 
these components is a maximal UV forest. Then we construct in a similar way the 
UV forests for other L! — 1 enumerations of the set of lines. After this we leave only 
distinct maximal UV forests. 



4 



Way 2. Since we consider a IVI graph we include it into any maximal forest. Let 
us delete a line from it. The resulting graph is decomposed as the union of its IVI 
components which we include into the maximal UV forest. Then we continue this 
process by deleting a line from some IVI component which is not a single line, etc. 

In the sector corresponding to a given maximal UV forest /, the function U takes 
the form 

76/ 

where Pf is a non-negative polynomial and the product is over elements of the given 
maximal UV forest /. We will call such a factorization proper. 

This factorization formula is proved by constructing an appropriate tree. One 
uses the relation 

= J]tM7)+'=(7nT)-cW^ (9) 

l^T 76/ 

where T is a tree or a 2-tree so that the factorization reduces to constructing a tree 
that provides the minimal value of the non-negative quantity 0(7 H T) — 0(7). Let Tq 
be the tree composed of all trivial elements of the given maximal UV-forest F. In 
other words, this tree can be constructed as follows. One uses an order of lines which 
was used within Way 1 for the construction of the given maximal UV forest / and 
includes the given line in the tree if a loop is not generated. One can obseve that this 
tree Tq provides the zero value of 0(7 fl Tq) — 0(7) for all the elements of the given 
maximal forest. 

To analyze convergence of the integral ([1]) large values of ai (in particular, to 
reveal infrared (IR) divergences) one has to take into account the exponent as well. 
A possible way is to separate the integration over every ai into (0, 1) and (1, 00) and 
then to deal with each of these 2^ regions separately — see, e.g. [TOl [11^. This can 
be enough for a general analysis but cannot be good from the practical point of view 
because the number of the resulting sectors will be too large. A more reasonable 
approach is to turn [5J to an integral with a compact integration domain, where both 
UV and IR divergences are somehow mixed up and manifest themselves as divergences 
at small values of parameters of integration. We will do this in the next section. 

3 Strategies to reproduce Hepp and Speer sectors 

In [131 , the starting point is the alpha representation (jl]) where primary sectors are 
introduced. The set of primary sectors corresponds to the different choices of a line 
in the given graph. At this step, one chooses a line / = 1, ...,L and defines a sector 
A; by ctj < a/, i ^ I and the sector variables by ai = Uai, i ^ I. The integration 
over ai is then taken due to the delta function whose argument is supposed to be the 
sum of all the variables minus one. 
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One can also start directly from ([T]) and introduce primary sectors < a^, i ^ 
I there. For example, in the case of I = L, using the homogeneity properties of 
the functions in the representation and explicitly integrating over a^;, we obtain the 
contribution of Ax, as 



{tTT''/^)''T{a-hd/2) 



'0 ^0 

^a-{h+l)d/2 



lL-1 

n 



a 



ai~l 



L-l 



m, 



-| a-hd/2 



1 ; 



where 



U = U{ai, ai-i, 1) , V = V(ai, . . . , a^^i, 1' 



(10) 



Without loss of generality let us consider only this primary sector. 

Let us remind that, for the case of non-zero masses, a general analysis of the 
factorization is not known even for Euclidean external momenta. Let us therefore 
turn to the pure massless case, as in [5J. Let us describe how sectors of Speer type 
can be introduced in such a way that the whole integrand of (JTOl) has a proper 
factorization. As we could see in the previous section, the use of /-sectors provides 
a proper factorization ([8]) of the function U so that the factor ^a^-('^+i)<^/2 f |io|) is 
properly factorized. However, these sectors generally do not provide a factorization 
of the second non-trivial factor. This can be seen using our example of Fig. [H We are 
going to use smaller sectors which are in fact obtained from the /-sectors generated 
by the graph F by a further decomposition. 

Let F°° be the graph obtained from F by adding a new vertex and connecting 
it with all the external n + 1 vertices by additional lines. These lines are only auxiliary 
and no propagators correspond to them. When writing down the function U for F°°, 
let us include, by definition, these additional lines into any tree. Then in the case of 
two external vertices (i.e. for n = 1) we have 

Vr = L{r°°Q^ 

where q is the only external momentum. 

Let us define sectors in a way similar to the previous section but using, instead of 
IVI subgraphs, another set of subgraphs which we call s -irreducible. If a subgraph 7 
does not have all the external vertices in the same connectivity component and if it 
is IVI let us call it s-irreducible as well. If a subgraph 7 has all the external vertices 
in the same connectivity component let us call it s-irreducible if the graph 7°° is IVI. 
Then the maximal forests consisting of s-irreducible subgraphs can be constructed 
again by Way 1 or Way 2. 

We define sectors (we name them Speer sectors) in a way similar to the sectors 
discussed in the previous section. We introduce sector variables by the same formula 
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([6]) as above. The factorization of the function V follows from its definition ([3]) and the 
auxiliary relation iQ. The 2-tree that provides a minimal value of the non- negative 
quantity 0(7 (IT) — 0(7) can be constructed by a procedure similar to the procedure 
used for the function U : one considers the lines in the order used for the construction 
of the given /-forest by Way 1 and includes the given line into the 2-tree if a loop is 
not generated and if this is not the line whose inclusion would connect all the external 
vertices. 

By construction, for such a 2-tree Tq, we obtain 0(7 fl Tq) — 0(7) = 6'(7) where 
^(7) = 1 if the external vertices are connected in 7 and ^^(7) = otherwise. Hence 
we obtain a proper factorization 

where Py is a non-negative polynomial. 

Obviously, the Speer sectors can be obtained from those associated with the graph 
r by a further decomposition, so that the factorization of the function U in the 
corresponding variables also holds and has the form similar to (jH]) with the same 
exponents. 

Let us turn to the modern strategies of sector decompositions. After introducing 
primary sectors, one obtains the contribution fllOp and other L — 1 contributions of 
the same type. At each step, one chooses a subset of the indices z/ = {ii, . . . 
and an index v from this subset and defines a sector < a^^, i ^ ir and the sector 
variables by = tjaj^, i ^ ir To formulate a strategy of introducing iterative sectors 
one needs to fix rules for determining subsets u at every recursive step 0. 

The first know sector decomposition strategy is described in [11]; three strategies 
guaranteed to terminate {A, B and C) and one strategy not guaranteed to terminate 
(X) are presented in [19]; they all are based on analyzing the functions U and V and 
choosing a subset of indices depending on their properties. Strategy S [20] is a bit 
different and is based on analyzing the polytopes of weights and their lowest faces. 

With the Hepp sectors, the situation is obvious: they are reproduced when we 
consider maximal subsets of lines at each step, i.e. with one line less than before this. 

To reproduce the choice of Speer sectors within a sector decomposition let us 
remind the Way 2 to construct sectors. One has to consider only subsets of the 
indices v = {ii, . . . ,ik} that correspond to s-irreducible subgraphs. 

We compared the number of sectors for numerous examples and discovered a 
curious phenomenon: the set of the sectors within Strategy S and the set of the Speer 
sectors is always the same. In particular, this is observed in the two examples of 
massless vertex diagrams shown in Fig. [2] at Euclidean external momenta, and in 
rather non-trivial examples of three four-loop propagators diagrams of Fig. [3] which 
are the most complicated master integrals among all four-loop massless propagators 

^Some strategies use more general sectors by comparing the integrations variables in different 
powers. 
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v2 \ v3 
Figure 2: Vertex off-shell diagrams 




m61 m62 m63 



Figure 3: Most complicated four-loop propagator diagrams 



integral^. The results of this analysis are shown in Table 1 where the number of 



diagram 


S 


X 


v2 


102 


102 


v3 


2160 


2251 


m61 


26208 


32620 


m62 


26304 


27540 


m63 


27336 


failed 



the sectors is shown. The first column stands for Strategy S and Speer sectors; 
the resulting number of the sectors is the same. For comparison, we included the 
second column for Strategy X [H] which has been proved to be very effective in a 
number of complicated calculations. Here 'failed' means that a factorization by sector 
decomposition has not been achieved for a reasonable amount of time (at least not 
for one day.) 

4 Conclusion 

In fact, our motivation to recall Speer sectors was to suggest to use them within 
FIESTA since they are optimal for Feynman integrals at all Euclidean external mo- 

^ Analytical results in expansion in e up to £° for these diagrams will be published soon [22] . 
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menta. However it turned out that these sectors are reproduced, in rather nontrivial 
examples, within Strategy S. Therefore we can conclude that Strategy S has chances 
to be an optimal universal strategy which always terminates and provides a proper 
factorization. 

Let us remind that Speer sectors have a certain physical meaning: the integration 
over the sector variable is responsible for UV divergences and, if 7 contains all the 
external vertices, for off-shell IR divergences. If on-shell and coUinear divergences are 
present, Speer sectors are no longer applicable. Then one can use sector decompo- 
sitions, e.g. within Strategy S, and this strategy could help to reveal the physical 
meaning of the sectors. To do this one might start with analyzing simple typical 
diagrams with on-shell or/and coUinear divergences. 
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